Relation between positional specific heat and static relaxation length: Application to 

supercooled liquids 



S. Davatolhagh 

Department of Physics, College of Sciences, Shiraz University, Shiraz 71454, Iran 

(Dated: February 2, 2008) 

A general identification of the positional specific heat as the thermodynamic response function 
associated with the static relaxation length is proposed, and a phenomenological description for the 
thermal dependence of the static relaxation length in supercooled liquids is presented. Accordingly, 
through a phenomenological determination of positional specific heat of supercooled liquids, we 
arrive at the thermal variation of the static relaxation length £, which is found to vary in accordance 
with £ ~ (T — To) - " in the quasi-equilibrium supercooled temperature regime, where To is the 
Vogel-Fulcher temperature and exponent v equals unity. This result to a certain degree agrees with 
that obtained from mean field theory of random-first-order transition, which suggests a power law 
temperature variation for £ with an apparent divergence at To. However, the phenomenological 
exponent v = 1, is higher than the corresponding mean field estimate (becoming exact in infinite 
dimensions), and in perfect agreement with the relaxation length exponent as obtained from the 
numerical simulations of the same models of structural glass in three spatial dimensions. 

PACS numbers: 64.70.Pf, 76.60.Es 



The deepest and most interesting unsolved problem in the theory of solids is probably that of the theory and the 
nature of glass and glass transition |]|. One of the most striking features of the typical supercooled liquid is that 
its relaxation time or viscosity r\ changes by several decades on changing the temperature by a few tens of degrees. 
All available data for viscosity fall between Arrhenius and highly non-Arrhenius extremes, designated 'strong' and 
'fragile', respectively Q. The latter is characterized by a highly temperature-dependent effective energy barrier 
against the viscous flow; i.e., a temperature-dependent energy barrier E e g appears in r\ = 770 exp(f3E e g), where tjq is a 
temperature-independent but species-dependent parameter of the orders of 1CP 2 -10 3 poise, f3 — l/fc^T, and ks is the 
Boltzmann constant. The temperature variation of the viscosity for fragile supercooled liquids is described accurately 
over a wide range of temperatures by the Vogel-Fulcher empirical equation jj] : 



The apparent divergence temperature T appearing in Eq. Q is called the Vogel-Fulcher temperature and is often 
found to be very close to the Kauzmann temperature Tk [J] , where a configurational entropy of the liquid extrapolates 
to zero This well known experimental fact, T « Tr is indeed a matter of considerable interest as it seems to 
suggest that the ideal glass transition temperature observed dynamically, and thermodynamically, must have common 
physical origin @. It should be also pointed out that in addition to non-Arrhenius variation of the viscosity with 
the temperature in the supercooled temperature regime, fragile liquids are also characterized by a distinct jump in 
the second-order thermodynamic functions such as the specific heat C p , the isothermal compressibility kt, and the 
thermal expansion coefficient a at the laboratory or calorimetric glass temperature T g , where T g > Tq. 

Furthermore, fragile supercooled liquids are also distinguished by a highly non-exponential relaxation response as 
they approach equilibrium when perturbed. Often the Kohlrausch- Williams- Watts (KWW) 0, or stretched expo- 
nential function is used to characterize the relaxation response of fragile liquids 4>(t) = exp[— (t/r)P], where (3 is the 
non-exponentiality parameter such that < j3 < 1, r is a relaxation time, and both are found to be temperature- 
dependent. The above relaxation response typical of fragile liquids is well explained in terms of the existence of 
dynamically heterogeneous regions in the supercooled liquid such that the relaxation in a given region is exponential 
but the average relaxation time r varies with a broad distribution among regions . This dynamical heterogeneity of 
the supercooled liquids has been further confirmed by the recent numerical |9( , and experimental research |lfj| , which 
confirm the presence of such regions with characteristic lengths spanning 100s of molecules. The above numerical 
and experimental confirmations of the dynamical heterogeneity of the supercooled glass forming liquids, lends further 
support to the notion of cooperatively rearranging regions (CRRs) in a supercooled liquid . 

There are several indications that the viscosity, and the various other structural relaxation times of a supercooled 
liquid must be correlated with the average size of a CRR, which is a concept dating back to the considerations 
of cooperative relaxation by the Adam and Gibbs [TJ- In their approach the increase of the effective potential 
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energy barrier E c g = zAfi, with Afi being largely the potential energy barrier against rearranging a single molecule 
in a cluster found by z molecules, is due to an increase in the cluster size z as the temperature is lowered. This 
cooperativity concept requires a characteristic static length £ characterizing the average linear size of a CRR. Evidently, 
an incremental increase in £ of the order of a few nanometers is magnified in an exponentially large (macroscopic) 
relaxation time as a consequence of which the supercooled liquid falls out of equilibrium on experimental time scales, 
hence, making any underlying static thermodynamic transition unreachable under the laboratory conditions. From a 
theoretical point of view, the mean field theory of random-first-order transition (or discontinuous spin glasses) that 
exhibits qualitative features in tandem with the structural glass phenomenology, suggests an approximate power law 
temperature dependence for the static relaxation length such that £ ~ (T — To) -1 ', where the apparent di verg ence 
temperature T is the Vogel-Fulcher temperature, and the exponent v = 2/d or v = 2/3 in d = 3 dimensions [Tall4||. 
There has been an empirical attempt to investi gate the temperature variation of £ the results of which are more or less 
consistent with the above proposed power-law Results reported for fragile liquid o-terphenyl are v — 0.69 ± 0.06 
with To = 203 ± 6 K. But the relevance to the static relaxation length of the experimental procedure adopted in 
Ref. [13, and the various interpretations of the experimental data have been called into question by the subsequent 
experimental investigations |l6j |. More recently, numerical simulations in three dimensions of the microscopic models 
that exhibit random-first-order transition in the mean field such as, the p-spin glasses |l7j . and the frustrated Ising 
lattice gas model |18| . are found to be more in favour of a static relaxation length exponent v = 1. 

In an attempt to clarify some of the discrepancies concerning the precise nature of the thermal dependence of 
the static relaxation length £ of the fragile supercooled liquids, as alluded to in the above discussion, we adopt a 
phcnomcnological approach to obtain the temperature variation of £. The main ingredients in this semi-empirical 
approach are: (i) the temperature-dependent potential energy barrier E c g against the viscous flow that is embodied 
in the Vogel-Fulcher equation for the viscosity, (ii) The thermodynamic response function bond susceptibility \b as 
applied here to the case of liquids, which is to be regarded as the response function measure of tendency for bond 
ordering or correlated relaxation of bonds into their low-lying energy states, brought about by the rearrangement 
of a CRR. Bond susceptibility was introduced earlier in the context of a diluted-bond model system relevant to the 
problem of glass transition, the thermodynamic properties of which were investigated by means of the Monte Carlo 
simulation 19] . Here, the very concept underlying bond susceptibility, i.e., correlated ordering or relaxation of bonds 
where intermolecular bonds are treated as distinct objects possessing internal degrees of freedom or energy states, 
is generalized and applied to the case of laboratory liquids. This approach paves the way for identification of the 
interaction or positional specific heat Cj as the thermodynamic response function associated with the characteristic 
length of relaxation £. Subsequently, a semi-empirical determination of the positional specific heat Cj for the general 
class of fragile liquids, is used to arrive at the thermal variation of static relaxation length £ that, by definition, gives 
the average linear size of a CRR in the liquid. 

The rest of this paper is organized as follows. In Sec. [D] we describe the relevant conceptual and theoretical back- 
ground concerning the various competing ordering processes in a liquid, and give a definition for the bond susceptibility 
Xb as a response function measure of tendency for correlated bond ordering. Sec. IIIII contains the derivations of the 
various relationships among thermodynamic and correlation functions relevant to the present discussion. As it will 
become evident in Sec. IIII Al there exists a relationship of the form Xb = T ■ Cj relating the bond susceptibility 
of a liquid in canonical ensemble with the interaction or positional part of the specific heat. Furthermore, in Sec. 
IIII Bl bond susceptibility is shown to be intimately related to static relaxation length such that essentially \b ~ /3£ 2 - 
These results essentially point at association of positional specific heat Cj as a thermodynamic response function with 
characteristic length of relaxation £, which is a novel concept brought to light in section ITlI CI This association of C 
and £ is then applied in Sec. HVI to develop a phenomenological description for the temperature variation of the static 
relaxation length in fragile supercooled liquids. Concluding remarks and a summary of the main results are presented 
in Sec. [V| 

II. RELEVANT BACKGROUND 

A. Two-order-parameter description of liquids 

Attempts have been made to incorporate frustration arising from the local ordering of bonds in a supercooled liquid 
through the introduction of a local order parameter characterizing the energetically favoured local arrangements of 
the liquid molecules, which are not consistent with the crystallographic symmetry favoured by the density ordering 
or crystallization. In this two-order-parameter description of the liquids due to H. Tanaka |20l |21| . the frustration 
arises from competition between density ordering and local bond ordering, explaining why some molecules crystallize 
easily without vitrification, while others easily form glasses without crystallization. The effect of density ordering is to 
maximize the density of molecules favouring a close-packed crystallographic symmetry, while local bond ordering tends 
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to improve the quality of bonds by reducing the bond energies at the local level. This model therefore emphasizes 
that introduction is necessary of a bond order parameter, in addition to the density p(r), in order to have a complete 
thermodynamic description of the liquid state, and, in particular, of the supercooled glass forming liquids. The 
energetically favoured local structures, such as, e.g., the icosahedral arrangements favoured by the spherical molecules 
|22| . are taken to be randomly distributed in a sea of normal liquid. It is further argued that the local structures 
with finite, but long life times, act as impurities and produce the effects of fluctuating interactions and symmetry- 
breaking random fields against density ordering in a liquid, in much the same way as magnetic impurities frustrate 
magnetic ordering in a spin glass system 23] . In this two-order-parameter description of the liquids the 'bond order 
parameter' S(r) is taken to be defined by the local concentration of the energetically favoured structures, and the 
average concentration of local structures S is estimated to be given by S ~ gs / g p exp[f3(E p — E s )], where E t and 
gi are the energy level and the number of degenerate states of the i-type structure. (i=p corresponds to the normal 
liquid while i=S to the energetically favoured local structures.) Thus, active bond concentration S(r) is taken to have 
a frustrating influence on crystallization at the local level, and each molecule intrinsically has the cause of disorder 
and random fields against the density ordering. 



Bond susceptibility is defined as a response function measure of tendency for bond ordering or correlated relaxation 
of bonds into their low-lying energy states, brought about by the rearrangement of a molecular group/CRR. Bond 
susceptibility apart from normalization is defined by 



where, (Mb) denotes the thermal-averaged bond energy order parameter characterizing the configurational energy of 
the system (more of which in Sec. and the average field Hb that is referred to as the bond ordering field is a 
self-generated molecular field favouring the local ordering of bonds and against the density ordering or crystallization 
1 1 Of . The above physical quantities are introduced in order to be consistent with the above two-order-parameter 
description of the liquids that recognizes two competing ordering processes in a liquid, namely, global density ordering 
that results in crystallization, and local bond ordering that is responsible for glass transition. Bond susceptibility is 
further expressed in terms of the equilibrium fluctuations of the bond energy order parameter 



where (nib) = (Mb)/N is the normalized bond energy order parameter characterizing the configurational energy, 
8mb = (nib — (mb)) is the corresponding fluctuation, ./V is the system size, and angular braces denote the usual 
thermal average. Eqs. and for bond susceptibility can be readily derived from the thermodynamic relation 
dG = —SdT — (Mb) dHb, which gives the chan ge in free energy G(T,Hb) of a system undergoing bond ordering 
as opposed to density ordering or crystallization |l9j . Evidently, the bond ordering field Hb is the thermodynamic 
conjugate- field that couples to the bond energy order parameter (Mb) (which characterizes the configurational energy), 
and can be regarded as the average concentration of energetically favoured local structures S. 

III. RELATIONS AMONG THERMODYNAMIC AND CORRELATION FUNCTIONS 

In this section the very concept underlying bond susceptibility, i.e., bond ordering or correlated relaxation of bonds, 
is generalized and applied to the case of laboratory liquids where we treat intermolecular bonds as distinct objects 
possessing internal degrees of freedom or energy states. The line of reasoning presented culminates in identification 
of positional specific heat as the thermodynamic response function associated with the static relaxation length. 



The bond energy order parameter (ni b) is defined as a measure of the bond-order prevailing in a system, and 
characterizes the configurational energy |l9j. By definition, it assumes large values when intermolecular bonds are 
in their low-lying energy states as for a bond ordered low temperature phase such as the glass, and is negligible 
when bonds are distributed uniformly among all possible energy states that is indeed the case when the thermal 



B. Bond susceptibility 




(2) 



Xb = N/3(5m 2 b ) , 



(3) 



A. Bond susceptibility and positional specific heat 
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energy is far in excess of the typical intermolecular binding energy. As a result, the bond energy order parameter of 
a liquid in canonical (NVT) ensemble can be simply defined in terms of configurational energy of the liquid. With 
$(ti,T2, ■ • • ,Tn) denoting the potential energy function of a liquid composed of TV molecules, the bond energy order 
parameter for this system is defined by 

(m b ) = -($)/JV. (4) 

Eq. (@J) satisfies all that is required of bond energy order parameter. On substituting this expression into the 
fluctuation-dissipation equation for the liquid in question we have 

Xb =(3(S<i> 2 )/N (5) 

where, (5$ = (<f> — ($)). Another response function of interest and of immense relevance to the problem of the glass 
transition is the specific heat, where for a liquid in canonical ensemble may be expressed as a sum of two terms, 
a kinetic part Cfc, and an interaction or positional part Cj. The above distinction follows from the fact that the 
liquid Hamiltonian consists of two distinct parts: a kinetic energy part X}<=i Pi/^ m i covering the degrees of freedom 
associated with the molecular momenta, and a potential energy part $(ri,r2, . . . , r^v) containing the contributions 
to internal energy from interactions or positional degrees of freedom. The positional part of the specific heat is indeed 
the temperature rate of change of configurational energy: 

N dT ' (6) 

It can be readily shown that an expression for the positional specific heat in terms of the equilibrium fluctuations of 
the configurational energy is given by 

d/k B = /3 2 (<5$ 2 ) /N. (7) 

On comparing Eq. JSJ for the bond susceptibility of a liquid in canonical ensemble with Eq. Q for the positional part 
of the specific heat, we arrive at the following simple result: 

Xb = T ■ Ct . (8) 

It becomes evident from Eq. (JSJ) that the bond susceptibility of a liquid in canonical ensemble, characterizing the 
tendency for bond ordering or correlated relaxation of bonds, can be simply interpreted as the response function 
positional specific heat. 

It is noteworthy that the result expressed by Eq. JHJ is readily verifiable for certain lattice models such as the 
two-dimensional Ising model j2^|, and some impurity variants thereof |25j. where analytic solutions are available. In 
particular, the four-spin correlation functions w(r) = (oio^ovoy+i) — (o~xo~2)(cr r crr+i) that can be also interpreted as 
two-bond energy correlation functions, with &\02 characterizing the energy of a reference bond, while cr r cr r +i that of 
a bond in a different location in the system, when summed over r or all distinct pairs of bonds essentially produce the 
specific heat that is entirely interaction or positional for the aforesaid lattice models: J2 r w ( r ) — de/d(3—l + e 2 , where 
e = (o~xo~2) |25j . This apparent connection between bond susceptibility and two-bond energy correlation functions 
will be used extensively next to establish a quantitative relationship between bond susceptibility and static relaxation 
length. 



B. Bond susceptibility and static relaxation length 



In the context of the bond ordering picture, a CRR can be viewed as a correlated region of relaxing bonds. Thus, 
the correlation length of such a region of bonds can be regarded as the characteristic length of (cooperative) relaxation 
£. Following Eq. ©, the bond susceptibility of a liquid in NVT ensemble is expressed as 

Xb = £ [{Ml) {M h f} (9) 

where, V is the liquid volume, and Mb is the extensive bond energy parameter the thermal average of which is the 
bond energy order-parameter characterizing the configurational energy. For short-range molecular interactions, which 
is almost always the case, Mb can be expressed in terms of the volume integral of a microscopic bond energy density 
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mb{r) = X^iSi m u5{r — r,), where mj, characterizes the energy of the ith bond, and Nb denotes the total number of 
bonds in the system. Thus, for a d-dimensional system we can write 

M b = [ d d rm b (r) (10) 



where, the integral is evaluated over the liquid volume, and m b (r) characterizes the energy of an intermolecular bond 
situated at r. On substituting this expression into Eq. © and simplifying, we have 



Xb=P J d d r {(m b (r) m b (0)) - (ro&(r)) (m 6 (0))] . (11) 



In Eq. (|11[) . mf,(0) characterizes the energy of a reference bond and mb(r) that of a bond at a distance r from 
the reference one. The quantity in the square brackets of Eq. Ijll|l . is the two-bond energy correlation function 
Gf,(r) = ( 5m b {r) 5m b (0) } quantifying the spatial correlation of the fluctuations of the bond energy order parameter. 
For an isotropic system Gf,(r) = G b (r). Furthermore, if we take the spatial variation of Gb{r) to be of the form 

00 

&M~$§, (12) 

where £ is a characteristic length beyond which the correlation function rapidly vanishes; the bond susceptibility of a 
liquid and the characteristic length of relaxation are thus related by 

Xb ~ P t 44 = 13 (13) 



We must point out that in a previous work the functional form of Eq. I|12|) has been used for a similar spatial correlation 
function in the context of a defect theory of relaxation to successfully recover the generalized Vogel-Fulcher equation 
for the viscosity |2fj. There the relaxation is considered to be brought about by the movements of mobile defects 
whose spatial correlation is governed by Eq. (|T2)l . In the context of the bond ordering picture, Eq. (|T2|) is applied to 
the two-bond energy correlation function Gb(r) as the structural relaxation is now considered to be a consequence of 
correlated relaxation of bonds within a region whose average linear size gives the static relaxation length £. 

As a corollary of the result expressed by Eq. (|13fl , we note that a possible diverging bond susceptibility Xb at some 
finite temperature T* must necessarily imply a diverging static relaxation length £ at the same temperature. That is, 
if the variation with temperature of the bond susceptibility for a system of interest is found to be a power law of the 
form Xb ~ {T — J 1 *) -76 , then the temperature variation of the relaxation length must be also governed by a similar 
power law £ ~ (T — T*)~ v such that the exponents are related by the scaling relation 

76 - 2^ . (14) 

We must emphasize that the above result is consistent with a standard result of statistical mechanics, namely, 
a = (2 — r)')v, where a is the specific heat exponent, and rf is the power law decay exponent of the energy-energy 
correlation function [23. In our treatment leading to Eq. (|14H . however, we have taken rf = for a supercooled liquid 
system that is corroborated by the numerical simulations of various models of structural glass in three dimensions 
One can further identify 75 with a, as expected. As an application to disordered systems, the results 
established here will be used in Sec. II VI in a phenomenological description for the static relaxation length of the fragile 
liquids. 



C. Positional specific heat and static relaxation length 



It has now become evident that the bond susceptibility of a liquid can be expressed as Xb = T ■ Ci, where Ci is that 
part of the specific heat containing contributions from interactions or positional degrees of freedom. Furthermore, 
bond susceptibility Xb or indeed Ci are shown to be intimately related to the static relaxation length £ such that 
essentially Xb ~ Ci ~ £ 2 . With their thermal behaviours so closely correlated, we therefore propose the identification 
of positional specific heat Ci as the thermodynamic response function associated with the static relaxation length 
£. Hence, we must further emphasize the significance of the role played by the specific heat, and, in particular, the 
interaction or positional part of it, in the problem of the glass transition. Unfortunately not enough is known about 
the precise behavior of the specific heat near Tq and much less about the interaction part of it, from an experimental 
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point of view, as the supercooled liquid falls out of equilibrium on experimental time scales at kinetic glass temperature 
T g for the reasons pointed out in the Introduction. It is generally believed that the excess specific heat over crystal 
value AC P that is regarded to be due to a subset of positional degrees of freedom involving transitions between 
inherent structures (or metabasins) of the potential energy hypersurface, rises with the decreasing temperature in 
the supercooled temperature regime, and a hyperbolic form AC P oc 1/T has been assumed in conjunction with the 
Adam-Gibbs equation for the viscosity to recover the Vogel-Fulcher equation (2^]. However, a drastically different 
AC P has been also used to accurately account for the viscosity of silicate glasses [2^. Hence, in the forthcoming 
section where the preceding results will be applied to the case of fragile systems, a semi-empirical approach is adopted 
to estimate the positional specific heat of fragile supercooled liquids. 



IV. APPLICATION TO SUPERCOOLED LIQUIDS 



In this section we present a phenomenological description for the thermal dependence of the static relaxation length 
for the general class of fragile liquids, in an attempt to clarify some of the discrepancies that were referred to in 
the Introduction. As it turns out, the result obtained via this phenomenological approach to a certain degree agrees 
with that obtained from mean field theory of random- first-order transition, also referred to in Sec. [I] that suggests 
a power-law temperature variation for the static relaxation length in the quasi-equilibrium supercooled temperature 
regime with an apparent divergence at the Vogel-Fulcher temperature To 

The positional specific heat Ci is the temperature rate of change of the configurational energy (Eq. ©). For the 
case of fragile supercooled liquids where the E e s of Eq. (ffl) is to be largely interpreted as a temperature-dependent 
'potential energy barrier' against the viscous flow 12, 30], one can consider a relationship of the form 

(4>)(T) ~ -E eS (T) (15) 

where, (</>) = ($) /N is the normalized configurational energy. Eq. Ijl5|) expresses the average depth or minima of 
the potential energy hypersurface explored by the liquid at each temperature in terms of the height of the effective 
potential energy barrier against the viscous flow. It should be also mentioned that in the present analysis it is only 
the temperature rate of change of the above quantities that is of interest. The above equation simply indicates that 
the higher the energy barrier E e g, the lower are the minima and consequently the configurational energy (</>). Eq. I|15|) 
is also consistent with a potential-energy-landscape representation of a supercooled liquid according to which a liquid 
is progressively confined to the deeper minima of the potential energy hypersurface with the decreasing temperature 
whereby it becomes more viscous due to an increased potential energy barrier E c s against the viscous flow [sol ], 
or alternatively, a reduction in the configurational entropy |3lj . Although this interpretation of the dynamics of 
supercooled liquids in terms of increasing barrier heights with the decreasing temperature is not the only one found 
in the current literature [32j, it is the one we deem appropriate for the purposes of the present discussion. Hence, the 
interaction or positional part of the specific heat of the fragile liquids can be approximated by 

Ci = -^L . (16) 

Using the effective potential energy barrier implied by the Vogel-Fulcher equation QJ, i.e., E c s = AksT/ (T — To), we 
obtain 

Ak B TT 

Xb = T-C l = — — — j . (17) 
(T - T ) 2 

Clearly, Eq. I|17l) implies a power law temperature variation of the form Xb ~ Ci ~ (T — T))~ 2 for bond susceptibility 
as well as positional specific heat of the fragile supercooled liquids, with an exponent jb — 2. From the scaling relation 
Eq. (|14f> . the exponent v governing the thermal variation of the static relaxation length of the fragile supercooled 
liquids is thus given by v = 7(,/2 = 1. Hence, with the effective potential energy barrier embodied in the standard 
form of the Vogel-Fulcher equation, we obtain 

^(T-T))- 1 (18) 

where, To is the Vogel-Fulcher temperature, and exponent v is equal to unity. 

One may repeat the same analysis this time using the generalized form of the Vogel-Fulcher equation: 

" = vo cxp {w^w) ' (19) 
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with the identification, E e g = BksT/(T — To) 7 , where B and 7 are constant parameters. The special case of 7 = 1 
results in the standard form of the Vogel-Fulcher equation being recovered, however, different values for parameter 7 
can be also found in the literature [2a]. The following is the result obtained with this rather generalized form of E e g 
that includes 7 as an extra parameter: 

= T . r. . = 

(T - T Q p+ 



Xb = T.Ci= bT-{T-T Q )] . (20) 



In the limit T -> T , Eq. (gOJ gives 
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Xb(T -> To) ~ — _ ^ )7+1 . (21) 



Eq. H21(l implies, Xb ~ Cj ~ (T — To)^ 7 " 1 " 1 -*, and from Eq. (|14|) . the characteristic length exponent is now given by 
v = (l+7)/2. Evidently an accurate experimental measurement of the parameter 7 appearing in the generalized Vogel- 
Fulcher equation is essential for a precise determination of the exponent v through the phenomenological procedure 
presented here. 

A discussion of the observed difference between phenomenological [y=Y) and mean field value (i/=2/3) of the static 
relaxation length exponent of the fragile liquids appears to be in order at this stage. This difference can be attributed 
to the mean field nature of the theory of random-first-order transition that becomes exact in infinite dimensions, 
and is believed to have an associated upper critical dimension d u = 6, which is significantly higher than d = 3 
space dimensions of supercooled systems. One therefore expects that the mean field estimate of the relaxation length 
exponent becomes increasingly accurate as the number of space dimensions approaches the upper critical value of 
six. Thus, it appears that the mean field theory of random-first-order transition, which presents qualitative features 
analogous with the structural glass phenomenology, provides a lower bound estimate of the static relaxation length 
exponent of a supercooled liquid. Indeed, there are other instances where mean field theories return lower estimates 
of correlation length exponents. A prominent example is the mean field theory of the continuous phase transitions 
(including the Ising model) that by hyperscaling has an associated upper critical dimension of four, which returns 
one-half for the correlation length exponent that is again exceeded by the exact two-dimensional value (unity) , and 
reliable numerical estimates (0.63) for the corresponding three-dimensional system. Nevertheless, It is a matter of 
considerable interest that the phenomenological value of the static relaxation length exponent v = 1, is precisely 
the value obtained from the three-dimensional numerical simulations of the microscopic models that exhibit random- 
first-order transition in the mean field limit such as the p-spin glasses 01 > an d the frustrated Ising lattice gas model 

In this section, a simple model for configurational energy in terms of increasing barrier heights with the decreasing 
temperature has been used that, despite simplicity, is applicable to various types of fragile liquids with predominantly 
ionic, Van der Waals, hydrogen, or covalent bonding. It would be also interesting to look at certain specific models 
such as the Rosenfeld-Tarazona relation for the Lennard- Jones liquid |3^| , or indeed any model that can be used to 
distinguish between positional and kinetic contributions to the specific heat in the context of the present work |34j . 
That effort is deferred to another work to be presented in due course. 



V. SUMMARY 



For a liquid in canonical ensemble it is shown that the bond susceptibility and interaction or positional part of 
the specific heat are related by Xb = T • Cj. Furthermore, bond susceptibility and static relaxation length are 
found to vary as X b ~ PS, 2 - These relationships essentially point at the identification of positional specific heat as 
the thermodynamic response function associated with characteristic length of relaxation — a proposition that further 
emphasizes the significance of the role played by the positional specific heat in the problem of the glass transition, 
and is likely to find further applications in the theory of disordered systems, as applied here to the case of fragile 
supercooled liquids in a phenomenological description for the thermal dependence of static relaxation length in those 
systems. 

Through the phenomenological approach, the temperature variation of the characteristic length of relaxation for 
fragile supercooled liquids is determined to be governed by a power law £ ~ (T— To) - " implying an apparent divergence 
at the Vogel-Fulcher temperature T , which to a certain degree agrees with the corresponding result obtained from 
mean field theory of random-first-order transition in that the apparent divergence temperature is To m both these 
cases. However, the phenomenological exponent v is found to be unity that is higher than the corresponding mean 
field estimate, hence, favouring a stronger temperature dependence for the static relaxation length in the supercooled 
temperature regime. This difference can be attributed to the mean field nature of the theory of random-first-order 
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transition, as discussed above. It is indeed a matter of considerable interest that the phenomenological exponent v = 1 
is in perfect agreement with the corresponding value obtained from the three-dimensional numerical simulations of 
the same models of structural glass that exhibit random-first-order transition in the mean field limit. 
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